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Computer-Oriented  Calculus 

The  enclosed  booklet,  which  emphasizes  a "discovery"  approach  to  mastering  the 
fundamental  concepts  of  calculus  through  use  of  interactive  computing,  was  prepared  by 
Dr.  Fred  Bell  for  Project  Solo.  It  actually  consists  of  8 modules  in  sequence.  These 
modules  can  be  used  as  a supplement  to  any  standard  textbook  in  calculus.  Dr.  Bell  con- 
ducted some  interesting  evaluation  with  similar  material  at  the  beginning  college  level  at 
Cornell  University.  He  found  a considerable  improvement  in  understanding  of  concepts 
(e.g.  meaning  of  a derivative)  on  the  part  of  students  who  had  computer  usage  over  those 
who  didn't,  but  little  difference  in  mastery  of  fixed  techniques  (e.g.  d(x2)/ax=2x).  Com- 
ments on  these  units  by  teachers  and  students  are  solicited.  It  should  be  noted  that  some 
of  the  material  on  limits  would  also  be  useful  in  non-calculus  courses. 

Corn-Share  News 

Those  of  our  readers  who  are  interested  in  technical  matte x*s  might  write  to  Corn- 
Share,  Inc.  (Box  1588,  Ann  Arbor,  Michigan,  48106)  for  a copy  of  Yol.  4,  No.  2 of  the 
Corn-Share  News.  In  addition  to  the  article  on  NEWBAS1C,  we  would  like  to  call  your 
attention  to  the  announcement  of  the  Sigma  7 Commander  11  service  on  page  19.  An  op- 
portunity to  hear  about  selected  details  of  this  new  system  was  given  to  some  of  us  last 
year  at  Research  Park  in  Ann  Arbor;  our  impression  then  was  that  it  will  very  likely  be 
one  of  the  most  sophisticated  time-sharing  systems  to  come  out  of  the  R & D labs  in  some 
time.  It  is  not  too  far  out  of  line,  in  our  opinion,  to  consider  allowing  students  who  have 
upgraded  themselves  to  an  appropriate  level,  and  who  have  prepared  justifying  proposals. 


to  access  such  a system.  The  other  announcement  of  interest  is  the  60  character  per  sec 


ond  service.  Such  high  speed  would  be  of  interest  to  the  potential  administrative  users  of 
computing  in  a school.  . 

The  SIGCUE  Bulletin 

The  enclosed  article  presenting  the  case  for  extended  versions  of  BASIC  was  pre- 
pared at  the  invitation  of  Dr.  Karl  Zinn  for  use  in  INTERFACE,  which  is  the  official  pub- 
| lication  of  SIGCUE  (the  ACM  Special  Interest  Group  on  Computer  Uses  in  Education). 

There  will  be  other  invited  articles  discussing  the  pros  and  cons  of  this  topic  in  the  fall 
i issue.  We  would  recommend  this  publication  highly  to  all  educators.  Subscription  in- 
formation can  be  obtained  by  writing  to:  SIGCUE  bulletin,  109  East  Madison,  Ann  Arbor, 
Michigan,  48104. 

Newsletter  Renewal  Form;  Project  Solo  Film 


We  are  cleaning  up  our  newsletter  mailing  list,  and  request  that  you  complete  the 


enclosed  form  if  you  wish  to  receive  our  mailings  during  the  coming  school  year.  Also^a  29 


minute,  16mm.  color  film  (optical  sound)  on  Project  Solo  is  now  available  for  rental.  If 
you  are  interested  in  obtaining  this  film,  please  fill  in  the  appropriate  information  on  the 
renewal  form. 
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The  Case  for  Extending  BASIC  as  an  Educational 
Programming  Language 

Thomas  A.  Dwyer 
Department  of  Computer  Science 
University  of  Pittsburgh 

There  is  only  one  way  to  make  something  better,  and  that  is  to  change 
it.  It  is  important  however,  not  to  confuse  the  word  'change'  with  the  word 
'destroy'.  On  the  contrary,  intelligent  and  useful  change  is  characterized 
at  least  as  much  by  a concern  for  preserving  the  good  aspects  of  that  being 
changed,  as  it  is  with  revision.  It  is  often  preferable,  in  fact,  to  view 
change  as  an  augmentation  rather  than  a drastic  re- shuffling  of  the  original 
system. 

I believe  that  these  general  observations  provide  the  key  to  answering 
the  question  "whither  hence"  that  faces  the  educational  community  today 
when  considering  computer  languages  for  interactive  computing.  The  ques- 
tion is  an  important  one,  for  there  is  little  doubt  that  the  use  of 
such  computing,  especially  in  community  colleges  and  secondary  schools* 
is  about  to  rise  dramatically. 

Having  worked  quite  closely  with  secondary  school  administrators, 
teachers,  and  students,  it  has  become  Very  obvious  to  me  that  the  impor- 
tance of  being  able  to  communicate  computer-based  curricular  ideas  in  a 
single  "mother-tongue"  within  the  user  community  cannot  be  overestimated* 
While  one  could  argue  that  better  languages  could  have  been  chosen  for  this 
function  (my  own  choice  at  the  present  would  be  PIL-X,  a fourth  generation 
derivative  of  JOSS),  the  de-facto  educational  mother- tongue  of  BASIC  turns 
out  to  be  an  extremely  good  compromise* 
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I do  not  call  BASIC  a compromise  in  a derogatory  sense,  but  only  to 

indicate  that  there  is  a vast  spectrum  of  user  needs  even  within 

the  educational  community,  and  that  BASIC  was  deliberately  aimed  at  the 

'novice1  end  of  the  spectrum.  This  was  a very  wise  decision  indeed,  and 

it  accounts  in  great  part  for  the  imminent  growth  in  educational  computing 

referred  to  above.  It  has  been  our  experience*,  however,  that  the  novice 

til  til 

stage  has  a very  short  life  span  for  even  8 and  9 grade  students,  and 
that  being  able  to  capture  .the  real  potential  of  computing  for  the 
full  span  of  the  educational  eaqaerience  that  lies  ahead  of  such  youngsters 
is  critically  dependent  on  the  software  we  give  them.  Such  software  must 
incorporate  a growth  potential  that  matches  the  fantastic  expansion  in  sophis- 
tication that  takes  place  in  the  latter  years  of  any  person's  education. 

For  reasons  that  space  does  not  permit  listing,  I believe  that  it. 
is  also  ■ desirable  for  the  same  software  to  accomodate  other  educa- 
tional uses  (e.  g.  administrative). 

In  order  to  accomodate  such  growth  in  our  own  work,  we  have  gone 
the  route  of  developing  an  ejq>anded  (and  expandable)  version  of  BASIC 
called  NEW  BASIC. We  originally  experimented  with  the  idea  of  preserv- 
ing the  purity  of  an  "urtext"  version  of  BASIC  as  the  principle  language, 
placing  all  new  features  (e.g.  a CAI  scan  syntax)  in  a pre-processor,  call- 
ing the  combination  NEWBASIC/ CATALYST.  However,  at  the  user  level, 
learning  to  use  the  pre-processor  (we  insist  on  user  involvement)  proved 
to  be  far  more  distracting  than  learning  to  add  new  language  elements  to 
one's  repertoire.  In  practice,  most  of  the  new  features  were  incorporated 
as  functions,  which  means  that  the  syntax  of  BASIC  was  not  altered*  Where 
the  syntax  was  augmented,  or  where  things  like  multiple  statements  on  a 
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lino  were  permitted*  it  was  always  in  the  spirit  of  serving  the  novice.  This 
works  well  in  practice.  Even  the  youngest  students  understand  the  meaning 
of  something  like: 

10  IP  X>5  PRINT  10*X  ELSE  LET  X=5  GOTO  100 

To  think  that  the  beginner  needs  to  have  this  spelled  out  for  him  in  a "simple” 

form  such  as 

10  IF  X>5  THEN  40 
20  LET  X=5 
30  GOTO  100 
40  PRINT  10*X 

is  just  as  much  of  an  educational  mistake  as  forever  speaking  to  youngsters 
with  a one-syllable  word  vocabulary! 

t 

On  the  question  of  adding  new  functions,  1 would  certainly  hope  for  uni- 
formity in  name,  parameter  lists,  and  effect.  But  it  seems  very  defeating 
to  me  to  argue  that  nothing  should  be  done  until  a standards  committee  de- 
signs such  functions.  Committees  continue  to  be  inefficient  mechanisms;  I 
would  rather  see  them  used  to  arbitrate  differences.  For  example,  we  use 
a function  NUM(X)  which  produces  a truly  random  integer  :N  in  the  range 

J j 

1<=N<=X.  If  a standards  committee  could  argue  that  it  is  better  to  call  this 
function  RAN,  or  that  the  range  should  include  zero,  fine.  In  the  meantime, 

I know  that  we  are  getting  more  beginners  involved  in  more  interesting 
things  because  we  can  let  them  incorporate  a statement  like: 

100  PRINT  NUM<50)  FOR  1=1  TO  10 

as  part  of  a bigger  exploration,  instead  of  the  following  (which  io_  required 
on  most  other  systems  today): 

10  PRINT  "TYPE  A NUMBER" 

20  INPUT  N 
30  FOR  1=1  TO  N 
40  LET  X=RND (Z) 

50  NEXT  I 
100  FOR  1=1  TO  10 
110  PRINT  INT(50*RND(Z)+1) 

120  NEXT  I 


I 
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The  two  codings  in  this  last  example  provide  a revealing  contrast, 
analogous  to  that  between  statements  in  a higher  level  language  and  the 
corresponding  code  in  assembly  language.  The  reason  I would  always 
argue  for  the  higher  level  version,  is  that  it  induces  higher  level  use  of 
computing  by  the  man  whose  channel  of  communication  with  the  machine  is 
the  language. 

In  sum,  then^I  would  advocate  expanded  versions  of  BASIC^which  also 
allow  a simple  form  of  extensibility  by  liters  through  a true  function- defining 
feature.  The  "original”  form  of  BASIC  should  always  be  a subset  of  such  an 
expanded  version,  and  every  effort  should  be  made  to  standardize  the  form  of 
the  new  features.  The  disadvantage  of  such  an  expanded  language  is  that 
(1)  the  compiler  rapidly  becomes  complex,  and  (2)  vendors  will  tend  to  re- 
sist "compatibility"  eiforts.  The  answer  to  both  of  these  difficulties  is  a 
very  strong  stand  on  the  part  of  users,  and  especially  funding  agencies,  in 
insisting  on  the  primacy  of  the  educational  community  as  a benefactor  of 
the  best  efforts  of  computer  technologists. 


This  experience  is  based  on  two  years  of  work  with  "Project  Solo", 
an  experiment  concerned  with  the  use  of  interactive  computing  in 
secondary  schools.  A newsletter,  sample  curriculum  modules,  and 
other  documentation  on  the  project  are  available  from  the  author.  A 
film  documenting  the  first  year  of  the  project  is  also  available  for 
rental  (16mm.  color,  29  minutes,  optical  sound). 

2 

The  NEWBASIC  system  was  developed  in  cooperation  with  Com-Share, 
Inc.,  a commercial  time-sharing  company  based  in  Ann  Arbor,  Michi- 
gan. One  advantage  to  involving  a commercial  time-sharing  vendor  in 
an  educational  project  is  the  capability  for  immediate  transportability 
to  a large  number  of  cities  (55  in  this  case)  via  their  continental  net- 
work, as  well  as  abroad. 

3 

I also  advocate  expanding  user  interaction  capabilities,  but  this  is  a 
system-design  rather  than  language-design' consideration. . ;The  total 
NEWBASIC  system  incorporates  a number  of  such  capabilities. 
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Project  Solo  Newgletter  Subscription  Request  Form  - September  1971 

Please  continue  to  send  1 copy  of  the  Project  Solo  Newsletter  along  with 
selected  modules  to: 

NAME:  

INSTITUTION:  

ADDRESS:  

ZIP: 


Reader  Comments 

1,  Material  you  like  most,  or  would  want  more  of: 


i 

! 

I 

! 2,  Things  you  would  like  to  see  added  (or  deleted  ) to  the  mailings: 


3,  Other: 


Film  Request 

t 

The  undersigned  is  interested  in  renting  the  Project ‘Solo  film  (1  week), 
to  be  shown  to  audience (sj  totalling  people.  Please  send  rental  informa- 
tion to:  NAME:  

INSTITUTION:  

ADDRESS:  

ZIP:  
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PROJECT  SOLO 

Ve.paA.tme.nt  o£  ComputeA.  ScXewce 
UntveJUtty  o£  PX.tt6buA.gk  (15213) 
Module .6  0017,  0019,  0031,  0049,  7094, 

0104,  0106,  7122 . 
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SOME  HINTS  ON  THE  LIMITATIONS  OF  COMPUTERS  FOR  CALCULATIONS 

When  you  are  solving  some  of  the  problems  in  this  booklet 
you  may  encounter  difficulties  that  appear  to  be  caused  by 
inconsistencies  in  mathematical  methods  but  that  are  actually 
caused  by  the  manner  in  which  computers  store  and  process 
numbers.  Computers  can  not  handle  numbers  which  have  too 
many  significant  digits  or  which  cure  too  large  in  absolute 
value.  In  NEWBASIC  you  are  limited  to  numbers  having  ten 
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significant  digits  that  range  in  absolute  value  from  10 
to  10*76.  Problems  also  may  arise  when  computers  try  to 
divide  by  numbers  close  to  zero,  subtract  two  numbers  which 
cure  very  close  together,  or  add  a very  small  decimal  number 
to  a very  large  ;7hole  number.  One  of  the  jobs  which  many 
mathematicians,  computer  scientists,  and  engineers  work  on 
is  that  of  devising  ingenious  schemes  for  getting  around 
these  limitations  of  computers. 

The  following  examples  show  a few  of  the  kinds  of  math- 
ematical problems  that  will  be  solved  incorrectly  by  computers. 

Example  1;  The  problem  of  finding  2793  can  not  be 
solved  directly  on  a computer  because  279  3 is  too  large 
in  absolute  value  to  be  calculated  by  an  instruction  such 
as  >10  LET  Y=2.0f793.  This  statement  will  cause  the 
computer  to  print  out  X IN  EXP(X)  TOO  LARGE  LINE  10. 

Example  2;  8 + 10“ 13  is  8.000G000000001,  but  the 

instructions : 

>10  LET  X=8 . 0 
>20  LET  Y=1 . OE-13 
>30  LET  Z=X  + Y 
>40  PRINT  Z 
>50  END 

will  give  the  following  printout  for  Z. 

I 8 
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Example  3:  The  problem  here  is  to  find  x 

for  x=. 000000001.  /x+3  ~ ' -3 

The  program: 

>10  LET  X=. 000000001 

>20  LET  Y=X/(SQRT(X+9.0)-3.0) 

>30  PRINT  Y 
>40  END 

gives  6.247225158  for  the  value  of  Y and  this  is  not  the 
correct  solution  to  the  problem. 

However/  we  can  use  the  computer  to  solve  this  problem 
if  first  we  rationalize  the  denominator. 

x = /x+9  +3 

/x+9  *3 

The  program: 

>10  LET  X=. 000000001 
>20  LET  Y=SQRT (X+9 . 0)  +3.0 
>30  PRINT  Y 
>40  END 

gives  the  correct  answer  of  6 as  output. 
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LIMITS  OF  SEQUENCES,  LIMITS  OF  SUMS  OF  SEQUENCES,  LIMITS 
OF  A FUNCTION 


The  speed  and  accuracy  of  computers  does  make  them 
particularly  usefull  for  exploring  limiting  processes 
in  mathematics.  For  example  the  sequence  defined  by 
f{n)=ni/n  (for  n»l,2,3,***)  is  made  up  of  terms  close 
to  1 for  large  values  of  n.  The  computation  of  a 
sufficient  number  of  terms  having. large  enough  indices 
to  enable  us  to  be  reasonably  certain  that  n ^ ^ 

is  quite  tedious  if  we  use  pencil  and  and  paper;  however, 
the  simple  program: 

>10  INPUT  N 

>20  LET  F=Nf (1.0/N) 

>30  PRINT  N,  F 
>40  GOTO  10 

permits  us  to  compute  enough  terms  of  the  sequence  to  dis- 
cover that  1 is  a good  guess  for  the  limit.  The  fol- 
lowing table  shows  that  this  sequence  converges  to  its 
limit  rather  slowly. 


N 

F 

N 

F 

N 

F 

N 

F 

1 

1.000 

6 

1.348 

20 

1.162  * 

200 

1.027 

2 

1.414 

7 

1.320 

30 

1.120 

300 

1.019 

- 3 

1.442 

8. 

1.297 

40 

1.097 

400 

1.015 

4 

1.414 

9 

1.277 

50 

1.081 

500 

1.013 

5 

1.380 

10 

1.259 

100 

1.047 

1000 

1.007 

Using  a computer  to  "make  a sequence  approach  its  limit" 
can  help  to  clarify  our  understanding  of  limits;  however, 
we  must  exercise  extreme  caution  in  drawing  conclusions. 
Since  we  can  never  look  at  all  of  the  terms,  any  conclusion 
about  the  limit  of  the  sequence  must  be  considered  at  best 
an  educated  guess  and  is  always  suspect  until  we  can  find 
some  way  to  prove  our  guess. 


1.2 


Problem  l!  Write  programs  and  use  them  to  help  you  to  make 
guesses  about  the  limits  of  these  sequences. 

'/nr 

2n+I0 


A. 

lim 

n+°° 

B. 

Ai® 

C. 

lim 

n-H» 

D. 

Ai® 

E. 

Ai® 

F. 

lim 

n+« 

G. 

lim 

n-K» 

(Be  careful,  for  large  values  of 
large  for  the  computer.) 


(What  happens  if  n=l?) 


n, 


211  is  too 


Problem  2: 


(This  function  has  a limit,  but  for  very 
large  n the  computer  gives  incorrect 
values  for  fl+ltj  n* 

For  each  real  number  p,  f(n)=5F  determines 

n the 


a sequence  for  n=l ,2,3,***.  For  some  values  of 
sequence  converges  and  for  others  it  deverges.  Write  a 
program  that  can  be  used  to  assist  us  in  finding 

for  any  value  of  p.  Use  your  program  to  find  the  values 
of  p for  which  this  sequence  will  converge. 

Problem  3:  Make  up  some  interesting  sequences  and  write 
programs  to  help  you  discover  their  limits. 


The  problem  of  finding  the  limit  of  the  sum  of  the 
terms  of  a sequence  can  also  be  approached  with  the 
assistance  of  a computer.  For  instance  the  following 
program  will  sum  the  first  n terms  of  the  sequence 
‘1,  1/2,  1/4,  1/8,  ***  for  any  integral  value  of  n. 

>10  LET  S=1 
>20  INPUT  N 
>30  FOR  1=1  TO  N 
>40  LET  F=1.0/ (2. Of I) 

>50  LET  S=S+F 
>60  PRINT  I,  S 
>70  NEXT  I 
>80  END 
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For  n=10  we'  have: 


I 


S 


I 


s 


1 

2 

3 

4 

5 


1.5 

1.75 

1.875 

1.9375 

1.96875 


6 

7 

8 
9 

10 


1.984375 

1.9921875 


1.99609375 
1 . 998046375 
1.999023437 


From  this  table  it  appears  that  the  limit  may  be  2. 

Problem  4:  Write  programs  and  use  them  to  estimate  the 
limits  of  the  sums  of  these  sequences : 


(A)  1,  1/3/  1/9/  1/27/*** 

(B)  (1/2)  2 , (1/3)  2 , (1/4) 


1111 

(C)  0!,  IT,  2T,  3T, ' ' ' . 

(D)  1,  -1/2,  1/3,  -1/4,  1/5 

(E)  1,  1/2,  1/3,  1/4,  ,1/5, 


We  can  build  upon  our  procedures  for  exploring  limits 
of  sequences  and  limits  of  sums  of  sequences  to  study  the 
limit  of  a function  defined  on  any  subset  of  the  real 


have  any  real  value  except  1/  we  must  decide  upon  a 
way  to  let  x approach  1 in  order  to  obtain  useful 
data  about  this  limit.  One  of  the  many  possibilities 
is  shown  next. 

This  program  with  N=5 


10  INPUT  N 
15  PR. 

20  FOR  1=1  TO  N 
30  LET  D=10+(-I) 

40  LET  Xl=l-D 
50  LET  X2=l+D 

60  LET  Yl=(Xl+3+Xlf2-Xl-l)/(Xl-l) 
70  LET  Y2=(X2+3+X2+2-X2-l)/(X2-l) 


1.4 


80  PR.  XI, Y1 

90  PR.  X2,Y2 

100  NEXT  I 

110  END 

is  output. 

0.9 

3.61 

1.1 

4.41 

0.99 

3.960100001 

1.01 

4.040099998 

0.999 

3.99600101 

1.001 

4. O'*  4000991 

0.9999 

3.999600113 

1.0001 

4.000399887 

0.99999 

3.99996071 

1.00001 

4.00003929 

The  output  indicates  (but  certainly  doesn't  prove)  that 
the  limit  may  be  4. 

You  can  see  that  the  programming  skills  needed  to 
write  programs  for  exploring  limits  of  sequences  and 
functions  are  minimal  and  even  people  with  little  ex- 
perience on  the  terminal  can  soon  be  exploring  some 
rather  sophisticated  mathematics  concepts. 


Problem  5:  Write  programs  and  use  them  to  make  an 


educated 

guess  about 

the  following  limits: 

lim 

(A) 

x+0 

(l+xJ 

lim 

sin (x) 

* 

(B) 

x+0 

X 

lim 

x*-l 

(C) 

x+1 

x-T~ 

(D) 

lim 

x+0 

X 

/x+4  -2 

(See  Example 

3) 

(E) 

lim 

x+0 

(1+x) 1/x 

(You  may  run  intd  a 

problem  here.) 

(F) 

Compare  Problem 

5,  part  E to  Problem 

1,  part  G. 
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APPROXIMATING  AREAS 

Our  problem  now  is  that  of  finding  approximations  to 
areas  bounded  by  the  graph  of  a function  and  an  interval 
of  the  x-axis.  For  example  suppose  we  want  to  find  the 
area  bounded  by  the  graph  of  y=xT  and  the  interval 
[1/  25}  of  the  x-axis.  This  area  is  illustrated  in  the 
following  figure. 


The  first  approach  we  will  use  to  approximate  this 
area  is  that  of  constructing  rectangles , the  sum  of  whose 
areas  is  an  approximation  to  the  area  in  question. 

The  next  set  of  diagrams  illustrates  ways  of  building 
rectangles  to  get  increasingly  more  accurate  approximations. 


U 


2.2 


x 

With  two  rectangles  the  approximation  is  12(1)**  + 

12(13)^  ~ 53.27.  With  four  rectangles  it  is  6(1)'5‘  + 

X -L  l 

6(7)2  + 6(13)2  + 6(19)?  = 69.66.  With  eight  rec- 

tangles our  approximation  is  76.39. 

You  can  observe  from  the  shaded  areas  in  the  figures 
that  the  difference  between  the  area  and  our  approximation 
(the  error)  decreases  as  we  increase  the  number  of  rec- 
tangles. The  procedure  of  doubling  the  number  of  rectan- 
gles was  continued  to  8192  rectangles  to  obtain  the  fol- 
lowing table. 


TJuriber  of 

Approximation 

Number  of 

Approximation 

rectangles 

*to  the  area 

rectangles 

to  the  area 

2 

55.266 

128 

82.291 

4 

69.661 

256 

82.479 

8 

76.388 

512 

82.573 

16 

79.594 

1024 

82.620 

3 2 

81.148 

2048 

82.643 

64 

31.912 

4096 

82.655 

8192 

8 2.661 

We  could  have  constructed  some  of  the  rectangles  to 
extend  above  the  graph  and  the  widths  of  their  bases 
could  have  been  varied  as  shown  below. 
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The  next  program  /AREA/  is  an  example  of  the  "rectan- 
gle method"  (Eulers  method)  for  obtaining  a set  of  approx- 
imations to  the  area  bounded  by  the  graph  of  any  specified 
polynomial  of  degree  three  or  less  for  any  interval  of  the 
x-axis . When  you  use  the  program  you  are  free  to  select 
the  polynomial  and  the  interval.  After  the  program  has 
calculated  an  approximation  to  the  area  bounded  by  ax3 4  + 
1 bx2  + cx  + d over  [r,s]  four  options  are  available  to  you: 

1.  Use  more  rectangles  to  get  a better  approx- 
imation. 

2.  Approximate  an  area  bounded  by  the  same 
function  over  a different  interval. 

3.  Approximate  the  area  bounded  by  a 
different  function. 

4.  Quit  and  do  something  else. 

When  using  the  program,  after  you  get  the  message 
NOW  TYPE  EITHER  34,  44,  54,  OR  1000 

1.  Type  54  for  option  1. 

2.  Type  44  for  option  2. 

3.  Type  34  for  option  3. 

4.  Type  1000  for  option  4. 

A sample  run  together  with  the  program  listing  is 
shown  on  the  next  two  pages. 

A sample  run  of  the  program  /AREA/. 

THIS  PROGRAM  COMPUTES  APPROXIMATIONS  TO  THE  AREA 
ROUNDED  RY  THE  X-AXIS  AND  THE  GRAPH  OF  ANY  POLYNOMIAL 
FUNCTION  OF  DEGREE  THREE  OR  LESS  OVER  AN  INTERVAL  OF 
THE  X-AXIS.  THE  METHOD  USED  TO  APPROXIMATE  THE  AREA 
IS  TO  BUILD  RECTANGLES  WHOSE  SUM  WILL  BE  AN  APPROXIMATION 
TO  THE  AREA  WE  ARE  TRYING  TO  FIND. 

THE  POLYNOMIAL  FUNCTION  IS  OF  THE  FORM 

A*X»3  + B*Xt2  + C*X  ♦ D 

YOU  CAN  GIVE  ANY  VALUES  (INCLUDING  D>  TO  A,B*C*  AND  D. 
0.K.*  ASSIGN  NUMERICAL  VALUES  TO  A,  R>  C>  AND  D. 

A * 
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/AREA/  (continued) 


SELECT  AN  INTERVAL  CR#SJ  0N  THE  X-AXIS 
BY  ASSIGNING  VALUES  T0  R AND  S* 

R = 

?D 
S = 

? 2 

H0W  MANY  RECTANGLES  D0  Y0U  WANT  T0  BUILD  0N  CR#S3? 

N 

?4 


THE  INTERVAL  IS  C D #2  J* 

THE  NUMBER  0F  RECTANGLES  IS  4 * 

THE  FUNCTI0N  IS 

D X t 3 ♦ 1 Xt2  + D X + D * 

THE  APPR0XIMATI0N  T0  THE  AREA  IS  *****  1*7S  ***** 


N0W  TYPE  EITHER  34#  44#  S4»  0R  1 0DD* 

? 54 

H0W  MANY  RECTANGLES  D0  Y0U  WANT  T0  BUILD  0N  C,R#S3? 
N * 

? 8 


THE  INTERVAL  IS  C D #2  3* 

THE  NUMBER  0F  RECTANGLES  IS  8 
THE  FUNCTI0N  IS 

D X 1 3 ♦ 1 Xt2  + D X + D * 

THE  APPR0XIMATI0N  T0  THE  AREA  IS  *****  2*1875  ***** 


N0W  TYPE  EITHER  34#  44#  54#  0R  1 0DD* 

? 54 

H0VJ  MANY  RECTANGLES  D0  Y0U  WANT  T0  BUILD  0N  CR#S3? 
N » 

? 64 


THE 

THE 

THE 

INTERVAL  IS  £ D 

NUMBER  0F  RECTANGLES 
FUNCTI0N  IS 

# 

IS 

2 

64 

3* 

♦ 

D 

Xt3  ♦ 1 Xt2 

+ 

D 

X + 

D 

THE 

APPROXIMATION  T0  THE 

AREA 

is 

***** 

2*604492188 

***** 

\ 


O 

ERJC 


N0W  TYPE  EITHER  34#  44#  54#  0R  1 0DD* 
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A listing  of  the  program  /AREA/. 

10  PR.  "THIS  PROGRAM  COMPUTES  APPROXIMATIONS  TO  THE  AREA" 

12  PR.  "POUNDED  BY  THE  X-AXIS  AND  THE  GRAPH  OF  ANY  POLYNOMIAL" 

14  PR.  "FUNCTION  OF  DEGREE  THREE  OR  LESS  OVER  AN  INTERVAL  OF" 

16  PR.  "THE  X-AXIS.  THE  METHOD  USED  TO  APPROXIMATE  THE  AREA" 

18  PR.  "IS  TO  BUILD  RECTANGLES  WHOSE  SUM  WILL  PE  AN  APPROXIMATION" 
20  PR.  "TO  THE  AREA  WE  ARE  TRYING  TO  FIND-" 

24  PR.  "THE  POLYNOMIAL  FUNCTION  IS  OF  THE  FORM" 

26  PR. 

28  PR.  "A*Xt3  + B*Xt2  + C*X  + D" 

30  PR. 

32  PR.  "YOU  CAN  GIVE  ANY  VALUES  < INCLUDING  0)  TO  A#R#C.  AND  D. " 

34  PR.  "O.K.#  ASSIGN  NUMERICAL  VALUES  TO  A,.  B,  C#  AND  D." 

36  PR.  "A  s " INPUT  A 

38  PR.  "R  = " INPUT  B 

40  PR.  "C  * " INPUT  C 

42  PR.  "D  = " INPUT  D 

44  PR.  "SELECT  AN  INTERVAL  CR#S3  ON  THE  X-AXIS" 

46  PR.  "BY  ASSIGNING  VALUES  TO  R AND  S." 

48  PR.  "R  * " INPUT  R 

50  PR.  "S  * " INPUT  S 

52  IF  R»=S  GOTO  500 

54  PR.  "HOW  MANY  RECTANGLES  DO  YOU  WANT  TO  BUILD  ON  CR.  S3?" 

56  PR.  "N  * " INPUT  N 

57  PR. 

58  PR- 

60  LET  P*INT<N> 

62  IF  C N< 1 ) OR  < N<>P)  OR  CN»500>  GOTO  600 

6 4 LET  T * 0 

66  LET  W s <S-R)/N 

70  FOR  X = R TO  <S-W>  STEP  W 

72  LET  T * T + ABS<  A*Xt  3+B*Xt2  + C*X  + 0) 

74  NEXT  X 
76  LET  T » T*W 

78  PR.  "THE  INTERVAL  IS  C"|RI ". "I  Si "3  « " 

80  PR.  "THE  NUMBER  OF  RECTANGLES  IS  "INI"." 

82  PR.  "THE  FUNCTION  IS" 

84  PR.  A| "Xt 3 + "!Bl"Xt2  + "ICI"X  + "|DI"*" 

86  PR. 

88  PR.  "THE  APPROXIMATION  TO  THE  AREA  IS  *****"! Tl "**+**" 

89  PR. 

90  PR. 

104  PR.  "NOW  TYPE  EITHER  34#  44.  54.  OR  1000." 

106  INPUT  K 

108  IF  <K<»34>  AND  <K<>44)  AND  <K<»1000)  AND  (K<>54)  GOTO  700 

1 10  IF  K=34  GOTO  34 

1 1 1 IF  K=44  GOTO  44 

1 12  IF  Ks54  GOTO  54 

113  IF  K* 1 000  GOTO  1000 

500  PR.  "R  MUST  BE  LESS  THAN  S*  TRY  AGAIN." 

505  GOTO  44 

600  PR.  "N  MU$T  BE  A POSITIVE  WHOLE  NUMBER  " 

605  PR.  "LESS  THAN  500.  TRY  AGAIN."  - 
61-0  GOTO  54 

700  PR.  "YOU  GOOFED.  TAKE  IT  FROM  THE  TOP." 

710  GOTO  104  - 

1000  END  lO 
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Problem  6s  use  the  program  /AREA/  to  find  good  approximations 
to  the  following  areas: 


A. 

X2 

over 

[0,  4] 

B. 

X2 

over 

t-4,  4] 

C. 

X3 

over 

1-2 t 0] 

D. 

X3 

over 

[0,  2] 

E. 

X3 

over 

[-2,  2] 

F. 

2x3 

-3x2+x 

-1  ov 

over  [-4,  3] 


Problem  7 : Write  a program  and  use  it  to  approximate  the 
area  bounded  by  sin(x)  over  the  intervals. 

A.  [-7 r,  0] 

B.  [0,  tt] 

C.  [-7T,  TT] 
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3.1 


NUMERICAL  INTEGRATION 


The  method  of  constructing  rectangles  to  approximate 
areas  is  one  of  several  strategies  which  can  be  used  in 
numerical  integration.  Remember  that  we  can  define  the 
f ^ lim  ^ 

integral  L f(x)dx  to  be  ?_^(Ax^)f (x^) , where  the 

Ax^*s  are  subintervals  of  the  x-axis  which  partition  the 

interval  [a#  b]  and  the  f(x^)*s  are  respective  points  in 

these  intervals.  To  obtain  a set  of  approximations  to 

n 

f(x)dx,  we  need  to  compute  ?al  (Ax^  f (xi)  where  n,  the 

number  of  subintervals  into  which  [a,  b]  is  divided,  is 
increased  for  each  approximation  and  the  Ax^'s  are  decreased, 
in  effect  we  are  constructing  rectangles  of  base  Ax^  with 
height  f (x^) , for  i=l,  2,  3,  4,***,  and  summing  the  areas 
of  all  the  rectangles.  Since  some  of  the  f(x^)'s  could  be 
negative,  which  occurs  over  subintervals  where  the  graph 
of  f (x)  dips  below  the  x-axis,  some  of  our  areas  (Ax^)f(x^) 
will  be  negative. 

The  next  program  /INTEGRATE/  is  an  illustration  of  the 
trapezoid  method  for  approximating  the  definite  integral 
of  some  function  f (x) . In  this  method  we  use  approximating 
trapezoids  rather  than  approximating  rectangles  as  is  shown 
in  the  diagram. 


An  approximation  to 


b 

f (x) dx  using  rectangles . 


3.2 


fb 

An  approximation  to  L f (x) dx  using  trapezoids , 


After  you  have  loaded  the  program  /INTEGRATE/  do  not 
type  >RUN,  You  must  alter  the  program  by  supplying  it 
with  the  function  whose  definite  integral  you  want  to 
approximate.  This  function  must  be  entered  as  line  100  of 
the  program.  For  example  if  you  want  to  integrate  the 
function  y=x2  + 2x  -1  over  some  interval,  you  must  type 
>100  LET  Y=X+2  + 2,0*X  -1*0, 

Now  the  program  is  complete  and  you  can  use  it  by  typing 
>RUN«  The  dependent  variable  must  be  called  Y and  the 
independent  variable  X;  so  line  100  should  always  have 
the  form 

>100  LET  Y— F (X) 

where  F(X)  is  the  function  of  X which  you  are  going  to  study* 

After  obtaining  an  approximation  to  f (x)dx  you  have 
three  choices  when  the  program  calls  for  the  change  code 
by  typing  the  message 

TYPE  THE  CHANGE  CODE* 
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Choice  1: 
Choice  2: 
Choice  3: 


Type  30  to  get  a better  approximation  by 
increasing  the  number  of  trapezoids. 

Type  18  to  integrate  the  same  function 
over  a different  interval. 

Type  300  to  stop  the  program. 


If  you  type  300  you  can  either  stop  or  try  the  program 

on  a different  function  Y=G (X)  by  typing 

>100  LET  Y=G(X) 

>RUN 


Problem  g:  Use  /INTEGRATE/  to  obtain  a set  of  approximations 

^ — "v- ” 

to  / e dx.  Hake  a guess  as  to  the  value  of  this  integral. 

Jo  ■ , 


Problem  9i  Use  integrate  to  approximate: 
A.  Jx2  logi  0 (x)dx 


B.  ^ (x)dx 

c-  Ji*  iog1 


(x)dx 


f 1 

Problem  10:  Use  /INTEGRATE/  to  approximate Jq  log^  (x)dx. 
Some  ingenuity  is  required  here  because  log^(O)  is  undefined. 

Problem  11:  Sketch  a graph  of  y=x3-9x.  Use  /INTEGRATE/  to 
approximate  the  following  definite  integrals. 


A.  j 

£43  (x3-9x)dx 

r 0 

B.  J 

L3  (x3-9x)dx 

r 3 

c. 

(x3-9x)dx 

Jo 

1 

L43  (x3-9x)dx 

Interpret  the  results  of  A,  B,  C,  D relative  to  your  graph. 
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A sample  run  of  the  program  /INTEGRATE/. 

>100  LET  Y = X t2  - 2.0+X 
>RUN 

THIS  PROGRAM  COMPUTES  APPROXIMATIONS  TO  THE 
DEFINITE  INTEGRAL  OF  THE  FUNCTION  WHICH  YOU 
SUPPLIED  ON  LINE  100#  OVER  AN  INTERVAL  CA#B3* 

NOW  TYPE  A VALUE  FOR  A. 

7 I 

TYPE  THE  VALUE  OF  B. 

?3 

HOW  MANY  SUBINTERVALS  00  YOU  WANT  CA#B3  DIVIDED  INT0? 
? 16 


THE  INTERVAL  IS  C 1 #3  3* 

THE  NUMBER  OF  APPROXIMATING  TRAPEZOIDS  IS  16 

THE  APPROXIMATION  IS  ******  0*671875  ****** 

TYPE  THE  CHANGE  CODE. 

? 30 

HOW  MANY  SUBINTERVALS  00  YOU  WANT  CA#B3  DIVIDED  INT0? 
7 100 


THE  INTERVAL  IS  C 1 #3  3* 

THE  NUMBER  OF  APPROXIMATING  TRAPEZOIDS  IS  100 

THE  APPROXIMATION  IS  ******  0.696799998  . ****** 

TYPE  THE  CHANGE  CODE* 

? 30 

HOW  MANY  SUBINTERVALS  D0  YOU  WANT  £A#B3  DIVIDED  INT0? 
7 400 


THE  INTERVAL  IS  C 1 # 3 3* 

THE  NUMBER  0F  APPROXIMATING  TRAPEZOIDS  IS  400 

THE  APPROXIMATION  IS  ******  0.674174994  ****** 

TYPE  THE  CHANGE  CODE. 

7 500 

TPS  CHANGE  CODE  MUST  BE  18# 30#  OR  300. 

TYPE  THE  CHANGE  CODE. 

7 30 

HOW  MANY  SUBINTERVALS  DO  YOU  WANT  CA#B3  DIVIDED  INTO? 
7 500 


THE  INTERVAL  IS  £ 1 #3  3* 

THE  NUMBER  OF  APPROXIMATING  TRAPEZOIDS  IS  500 


THE  APPROXIMATION  IS  ****** 
TYPE  THE  CHANGE  CODE. 

7 


0*672671988 


£ + ££  + £ 


300 


A listing  of  the  program  /INTEGRATE/ 


10  PR-  "THIS  PROGRAM  COMPUTES  APPROXIMATIONS  TO  THE" 

12  PR.  "DEFINITE  INTEGRAL  OF  THE  FUNCTION  WHICH  YOU" 

14  PR*  "SUPPLIED  ON  LINE  100#  OVER  AN  INTERVAL  t A # B j * " 

16  PR* 

18  PR*  "NOW  TYPE  A VALUE  FOR  A*"  INPUT  A 
20  PR*  "TYPE  THE  VALUE  OF  B*"  INPUT  B 
25  IF  A»*B  GOTO  200 

30  PR*  "HOW  MANY  SUBINTERVALS  DO  YOU  WANT  CA#BJ  DIVIDED  INTO 
35  INPUT  N 

40  LET  N*FIX<ABS<N>> 

45  LET  S*0 
50  LET  W=(B-A)/N 
55  FOR  X=A  TO  B STEP  W 
100  LET  Y = Xt2  - 2*X 
105  IF  X * B GOTO  125 
110  IF  X s A GOTO  125 
1 1 5 LET  S=S  + 2*Y 
120  GOTO  130 
125  LET  S = S ♦ Y 
1 30  NEXT  X 

135  LET  T « < <B-A>/<2*N> >*S 
1 40  PR* 

142  PR* 

144  PR*  "THE  INTERVAL  IS  C "l A*"# "I B* "3 * " 

146  PR*  "THE  NUMBER  OF  APPROXIMATING  TRAPEZOIDS  IS"*N 
148  PR. 

150  PR*  "THE  APPROXIMATION  IS  ******"*  T*  ••******" 

1 52  PR* 

1 SS  PR.  "TYPE  THE  CHANGE  CODE*"  INPUT  C 
1 60  IF  C*18  GOTO  18 

165  IF  C*30  GOTO  30 
170  IF  Cs300  GOTO  300 
175  GOTO  250 

200  PR*"A  MUST  BE  LESS  THAN  B*  TRY  AGAIN." 

205  GOTO  18 

250  PR.  "THE  CHANGE  CODE  MUST  BE  J8#30#  OR  300*" 

2 55  GOTO  155 
300  EN0 


4.1 

APPROXIMATING  VALUES  OF  DERIVATIVES 

The  derivative  of  f(x)at  x=a  can  be  defined  as 

lim  f (ath)  -f  (a) . Since  a is  a constant,  f(ath)-f  (a) 
h+0  h T 

is  a function  of  h;  consequently  any  method  developed  for 
approximating  the  limit  of  a function  can  be  used  to  ap- 
proximate this  limit; 

To  approximate  the  derivative  of  y=x2for  x=3  directly 

from  definition,  we  need  to  find  values  of  (3th)2- (3) 2 for 

!E 

values  of  h which  are  close  to  0.  The  program  shown  below 

produces  the  following  set  of  approximations  to  lim  (3th) 2- (3) 2 

h+0  h 

where  N=5.  Notice  that  in  this  program  h is  approaching  0 
from  the  positive  side. 


10  INPUT  N 
20  FOR  1=1  TO  N 
30  LET  H=l/(10+N) 

40  LET  D=( ( (3tH) *2) - (3*2) )/H 
50  PRINT  H,  D 
60  NEXT  I 
70  END 


h 

(3+h) 2-  (3) 2 

h 

.1 

6.099999999 

.01 

6.009999994 

.001 

6.000999885 

.0001 

6.000099238 

.00001 

5.999999121 

Problem  12s  . Write  a program  which  will  compute  values  of 

(3th)  2-  (3)  2 Jjor  negative  values  of  h.  Use  your  program  to 
— h 

estimate  lim  (3th)-2-  (3) 2 as  h approaches  0 from  the  ne- 
h+Q  h 

gative  side.  What  is  the  derivative  of  x2  when  x = 3? 


4.2 


The  problem  of  finding  (3+h) 2- (3) 2 for  values  of  h 

h 

which  are  close  to  0 illustrates  one  of  the  difficulties 
that  can  occur  when  we  use  a computer  to  perform  our  cal- 
culations. Suppose  h-. 00000000001,  then  the  computer  must 
find  (3+. 00000000001) 2-32.  However,  3+. 00000000001  = 
.00000000001 

3.00000000001,  which  has  12  significant  digits,  and  our 
computer  can't  handle  this  many  significant  digits,  con- 
sequently the  computer  will  take  3.00000000001  to  be 

3.0000000000  and  compute  (3.0000000000) 2-(3. 0000000000) 2 

1.0E-11 

which  is  0.  The  correct  answer  is  6.00000000001,  and  in 

fact  we  are  going  to  get  erroneous  data  when  |h|<.0001. 

We  can  avoid  this  problem  by  simplifying  (3+h) 2-  (3) 2 before 

h 

we  write  our  computer  program. 

(3+h)2- (3) 2 _ 9+6h+h2-9  _ g+h 
h — h ~ 

So,  lim  (3+h)  2-  (3)  2 __  lim  (6+h) , and  for  h close  to  0, 
h+0  Ii  h+0 

6+h  is  close  to  6. 

If  we  don't  make  this  simplification  the  computer  gives 

us  incorrect  information.  If  we  make  the  simplification  we 

don't  really  need  a computer  to  show  us  that  lim  (6+h)  = 6. 

h+0 

In  this  problem  the  computer  really  isn't  of  much  help. 


As  a general  rule  when  using  the  computer  to  find 

f (x+h)  -f  (x)  for  various  values  of  x,  you  should  avoid 
h 

letting  h be  closer  than  .00001  to  0 as  you  attempt  to 

estimate  lim  f (x+h)-f (x)  . 
h+0  ~E 


Problem  13:  Write  a program  which  can  be  used  to  find  esti- 
mates for  the  derivative  of  log  x for  any  value  of  x whicl 

io 

the  user  of  the  program  chooses.  You  will  need  to  consider 

finding  lim  log  (x+h) -log  (x) . 
h+0  — 5 — 
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4.3 


Remember  that  log10x  is  only  defined  for  positive  values  of  x. 
Test  your  program  by  estimating  the  derivatives  for  x=10,  100, 

1,  2,  and  5. 

The  next  program  called  /DERIVSIN/  is  one  which  can  be 

used  to  find  the  derivative  function  of  sin (x) . Our  strategy 

will  be  to  find  approximations  to  lim  sin(x+h)-sin(x)  for  manv 

h+0  h 

values  of  x,  obtaining  a set  of  pairs  of  numbers  {(x,  approxi- 
mation to  lim  sin  (x-frh) -sin (x)  vifor  these  values  of  x.  We  will 
h+0  h * 

then  plot  all  of  these  points  on  a rectangular  coordinate  system, 
sketch  the  graph  which  they  form,  and  try  to  guess  the  function 
which  has  this  graph.  For  each  value  of  x that  is  entered 
into  the  program,  it  will  give  as  output  jc,  sin(x)  and  values 
of  h and  sin (x-frh)-sin (x)  for  h-±.l,  ±.01,  ±.001,  and  ±.0001. 

n 

Problem  14:  Use  the  program  /DERIVSIN/  to  find  the  derivative 
function  of  sin (x) . Select  values  for  x until  you  have  enough 
approximations  to  the  derivative  of  sin(x)  to  sketch  reasonably 
good  graphs  of  both  sin (x)  and  D sin (x) . What  is  your  guess 
as  to  the  derivative  function  of  sin(x)? 

Problem  15:  write  a program  similar  to  /DERIVSIN/  which  for 
any  value  of  x will  give  as  output 

x,  ex 

and  h,  ex+h-ex 

h 

for  h=±.l,  ±.01,  ±.001,  and  ±.0001.' 

Use  this  program  to  sketch  graphs  of  y=ex  and  the  derivative 
function  of  ex. 

Problem  16:  write  a program  which  for  any  value  of  x will 
give  as  output  x,  cos(x)  and  h,  cos (x+h) -cos (x) 

for  h*±.l,  +.01,  +.001,  and  +.0001. 

Use  this  program  to  discover  derivative  function  of  cos  (x) . 
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A sample  run  of  the  program  /DERI VS IN/. 

PROGRAM  T0  ESTIMATE  THE  DERIVATIVE  0F 
SINCX)  F0R  ANY  VALUE  0F  X. 

WHAT  IS  X? 

?-1 


X « -1  SINCX)  = -0.841470985 

H C S INC X+H) -SINCX) 3/H 


0*1  0*581440752 

-0*1  0*497363753 

9.999999999E-03 
-9.999999999E-03 
9.999999999E-04 
-9.999999999E-04 
9* 999999998 E- 05 
-9.999999998E-05 


0*54450062 

0*536085982 

0*540722944 

0*53988149 

0*540344336 

0*540260335 


WHAT  IS  X? 
?0 


X * 0 SINCX)  = 0 


H C S INC X+H) -SINCX) 3/H 


0*1  0*998334166 

-0*1  0.998334166 

9.999999999E-03 
-9.999999999E-03 
9.999999999E-04 
-9.999999999E-04 
9.999999998E-05 
-9.999/99998E-05 


0*999983333 

0*999983333 

0*999999833 

0*999999833 

1 

1 


WHAT  IS  X? 

?1 

X s 1 SINCX)  « 0.841470985 


H C SINCX+H) -SINCX) 3/H 


0.1  0.497363753 

-0*1  0*581440752 

9.999999999E-03 
-9* 999999 999E-03 
9.999999999E-04 
-9.999999999E-04 
9.999999998E-05 
-9. 999999 998 E-05 


0*536085981 

0*544500621 

0.539881479 

0.540722951 

0*540260226 

0*540344372 
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WHAT  IS  X? 
?*-^SC* 


16 
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A listing  of  the  program  /DERI VS IN/. 

10  PR.  "PROGRAM  TO  ESTIMATE  THE  DERIVATIVE  OF" 
IP  PR.  "SIN<X>  FOR  ANY  VALUE  OF  X." 

14  PR. 

16  PR*  "WHAT  IS  X?"  INPUT  X 
18  PR* 

20  PR* 

30  LET  Z*SIN<X> 

40  PR*  ,#X  » "!XI"SIN<X>  a "IZ 
42  PR. 

SO  PR."  H CSIN(X+H)-SINCX) 3/H" 

S2  PR. 

60  FOR  N=1  TO  4 
70  LET  HI  a lO.Ot(-N) 

80  LET  H2  a -HI 

90  LET  Y1=<SIN<X+H1)-Z)/H1 

100  LET  Y2a<SlNCX+H2)-Z)/H2 

1 10  PR.  HU  Y1 

120  PR*  H2#  Y2 

130  NEXT  N 

140  GOTO  14 


4.6 


The  following  program  is  a generalization  of  the  pro- 
gram /DERI VS IN/  which  can  be  used  to  explore  the  derivative 
function  of  sin  (x) . When  using  this  program,  you  have  the 
option  of  selecting  any  function  whose  derivative  you  want 
to  study.  To  use  this  program  called  /DERIVTIV/,  first  load 
the  function  into  your  workspace.  Don't  try  to  run  this 
program  until  you  alter  it  by  providing  the  function  whose 
derivative  you  are  going  to  study.  You  can  do  this  by  typing 
>100  LET  Y=F(X) 

where  F(X)  is  any  function  you  want.  For  example  to  find 
approximations  to  values  of  the  derivative  function  of 
y=/x+4  you  must  type 

>100  LET  Y=SQRT(X+4.0) 

Of  course  you  must  remember  that  x must  be  greater  than  or 
equal  to  -4  or  SQRT (X+4.0)  will  be  imaginary. 

When  you  finish  using  the  program  to  explore  your  function 
hit  the  ESCAPE  key.  Now  you  can  either  LOGOUT  or  work  on  an- 
other function  y=G (X)  by  typing 

>100  LET  Y=G(X) 

>RUN 

where  G(X)  is  a specified  function  of  X. 

The  following  pages  contain  a set  of  sample  output  from 
/DERIVTIV/  and  a listing  of  the  program. 

Problem  17 : Use  /DERIVTIV/  to  obtain  enough  data  to  plot  an 
accurate  graph  of  y=A  and  the  derivative  function  of  y=i. 

Plot  these  graphs  and  try  to  guess  the  function  which  is  the 
derivative  of  y=i. 

Problem  18:  Use  /DERIVTIV/  to  assist  in  plotting  the  graphs 
of  y=|x|  and  its  derivative  function.  Try  to  'guess  the  der- 

t 

ivative  function  of  y=|x|. 

Problem  19:  Use  /DERIVTIV/  to  assist  in  plotting  the  graphs 
of  y=sin  (2x)  and  the  derivative  function  of  y=sin  (2x) . Can 
you  guess  what  this  derivative  function  is? 


3Q 
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A sample  run  of  the  program  /DERIVTIV/ . 


>100  LET  Y=EXP(X) 

>RUN 

THIS  PROGRAM  WILL  COMPUTE  A SET  0F  APPROXIMATIONS 
TO  THE  DERIVATIVE  OF  ANY  FUNCTION  F<X>  WHICH  YOU  GIVE 
IT  FOR  ANY  VALUE  YOU  ASSIGN  TO  THE  INDEPENDENT 
VARIABLE  X#  PROVIDED  THE  FUNCTION  HAS  A DERIVATIVE 
FOR  THAT  X. 

FOR  EACH  X YOU  WILL  GET  AS  OUTPUT  THE  VALUES  OF  X 
AND  F<X)  # AS  WELL  AS  A SET  OF  APPROXIMATIONS  TO 
F * (X)  FOR  H s .1.  -.1#  .01#  -.01#  .001#  -.001# 

. 0001 # AND  - .0001 • 

WHAT  IS  X? 

?0 


H APPROXIMATION  TO  F*<X) 

0. 1 1.051709181 

-0.1  0.951 62582 

9.999999999E-03  1 .005016708 

-9.999999999E-03  0*995016639 

1.000500161 
0 • 999500247 

1 .000049961 
0.999950571 


X = 0 Y * I 

WHAT  IS  X? 

? 1 

H APPROXIMATION  TO  F#<X) 

0*1  2 .858841955 

-0. 1 2.586787173 


9.999999999E-04 

9.999999999E-04 

9.999999998F-05 

9.999999998E-05 


9.999999999E-03 

9.999999999E-03 


2*73191865 

2*704735608 


9 • 9999999 99E-04 
-9.999999999E-04 


2.71964143 

2.716923176 


9.999999998E-05 

-9.999999998E-05 


2.718417673 

,2*71814628 

t * 

i-  - ■,  i 


X - 1 

WHAT  IS  X? 
?**ESCJ 


2.718281828 

31 


Y * 
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A listing  of  the  program  /DERIVTIV/. 


A PR.  "THIS  PROGRAM  WILL  COMPUTE  A SET  0F  APPROXIMATIONS" 
6PR.  "TO  THE  DERIVATIVE  0F  ANY  FUNCTI0N  F<X)  WHICH  Y0U  GIVE" 
8 PR.  "IT  F0R  ANY  VALUE  Y0U  ASSIGN  T0  THE  INDEPENDENT" 

10  PR.  "VARIABLE  X#  PR0VIDED  THE  FUNCTION  HAS  A DERIVATIVE" 

1  2 PR.  "F0R  THAT  X." 

1 A PR. 

16  PR.  "F0R  EACH  X Y0U  WILL  GET  AS  0UTPUT  THE  VALUES  0F  X" 

18  PR.  "AND  F<X  ) # AS  WELL  AS  A SET  0F  APPR0X I MATI 0NS  T0M 

20  PR.  "F  * < X ) F0R  H * .1#  -.1#  .01#  -.01#  *001#  -*0Q1#" 

22  PR.  ".0001#  AND  - .0001." 

2 4 PR. 

2 5 PR. 

26PR  • 

28  PR.  "WHAT  IS  X?"  INPUT  X 
30  PR. 

32  PR.  " H APPR0X I MAT 1 0N  T0  F’<X)" 

35  FOR  N® 1 T0  4 

40  LET  H=!0t<-N) 

45  LET  C * 0 

47  LET  1=0 

50  LET  A=X 

55  LET  X=X+M 

100  LET  Y=EXP<  X ) 

110  IF  I = ! 60T0  200 
120  LET  I =1  + 1 
I 30  LET  0 1 *Y 
I 40  LFT  X=A 

1 50  G0T0  100 
200  LET  02=Y 

210  LET  0=(0I-02)/H 
220  PR.  H#  0 
230  LET  C=C+I 
240  IF  C=2  GOTO  270 
250  LET  H = <-H) 

2 60  GOTO  47 
270  PR. 

300  NEXT  N 
320  PR. 

3 30  PR. 

340  PR.  "X  = "|Al  " Y = "102 

350  G0TO  25 


t 
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THE  FUNDAMENTAL  THEOREM  OF  INTEGRAL  CALCULUS 


5.1 


/. 


A modification  of  /INTEGRATE/  can  be  used  to  "numerically 
verify"  The  Fundamental  Theorem  of  Integral  Calculus;  i.e. , 
k f (x)dx=F (b) -F (a)  where  F (x)  is  a primitive  of  f (x) . 


For  example, 

sin(x)dx={^cos(Tr^  Tj-EcosCO)] 


t 


can  be  verified  by  calculating  a sequence  of  approximations 
sin(x)dx  and  comparing  the  suspected  limit  of  the 


to 


sequence  to  Ecos^  ^“[jcos  (0^|=1. 


The  following  program/  /FUNDTHEO/,  uses  the  trapezoid 
method  to  compute  approximations  to  f(x)dx  for  a given 
function  f (x)  over  a specified  interval  [a,  b]  for  n trap 
ezoidS/  and  compares  the  approximation  to  g(b)-g(a)  where 
g (x)  is  a primitive  of  f (x) « 


After  you  have  loaded  the  program , modify  lines  100  and 
200  as  follows: 

Line  100  should  contain  the  function  to  be  integrated/ 

Y=F (X)  in  the  form: 

>100  LET  Y=*F(X) 

Line  200  should  contain  any  primitive  Z=G(X)  of  Y=»F  (X) 
in  the  form: 

>200  LET  Z=G(X) 

For  example,  if  the  integral  cos (x)dx,  then  line 

100  should  be: 

>100  LET  Y=*C0S  (X) 
and  line  200  could  be: 

>200  LET  Z=SIN(X)+1 

since  SIN (X)+l  is  a primitive  of  COS  (X) . 

After  you  have  supplied  lines  100  and  200  you  can  execute 
the  program  by  typing 
>RUN 


o 
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5.2 


1. 


After  the  program  has  printed  an  approximation  to 
f(x)dx  and  G(B)-G(A)  it  will  print 


TYPE  THE  CODE 

To  use  a different  number  of  approximating  trapezoids  on 
the  interval  [a,  b]  type  30;  to  get  a new  interval  by  chang- 
ing the  values  of  a and  b,  type  20;  to  stop  the  program 
type  300. 


Problem  20:  Use  /FUNDTHEO/  to  approximate  ^-x“dx  an<j  to 
verify  that  jT^-x^dx=|-j  +2I  . 

° L r v*  2 

Problem  21;  Use  /FUNDTHEO/  to  approximate  L cos  (x)dx 


and  to  verify  that 


/. 


VA  2 

cos  (x)dx= 


x + sin (x) cos (x) 


2 * 2 
dx 


f 1 

Problem  22:  Use  /FUNDTHEO/  to  verify JQ 

Problem  23;  Use  /FUNDTHEO/  to  verify 

“4 


-3 


TT2 


= |arctan  (x)~|^ 


I2 


4 3x 
* eoxdx* 


e3x 

| -6 


r 


log  (x) dx= 
© 


x(logex)-x 


3 

2 


C.  sec“(x)dx=|  tan  (x) 


c=^tan  (x ) -2J  * 


D.  J ^ /1+cos  (x)  dx=j2/!T  (sin^J 
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A sample  run  of  the  program  /FUNDTHEO/ . 


5.3 


* RUN 

THIS  PROGRAM  COMPUTES  APPROXIMATIONS  TO  THE  DEFINITE 
INTEGRAL  OF  THE  FUNCTI0N  Y * FCX)  WHICH  YOU  SUPPLIED  0N 
LINE  100#  OVER  AN  INTERVAL  £a#BI|,  AND  C0MPARES  THE 
APPROXIMATION  T0  G<fe)-GCA'#  WHERE  Z = G<X>  IS  A 
PRIMITIVE  0F  FCX)#  WHICH  YOU  SUPPLIED  0N  LINE  200. 

WHAT  IS  A? 

? -3* I 41 529 
INPUT  B. 

?0 

HOW  MANY  SUBINTERVALS  0N  CA#BI? 

? 200 

INTERVAL  C -3*141529  #01.  200  TRAPEZOIDS 

APPROXIMATION  TO  INTEGRAL  IS  *******  -1.999958877  ******* 

G< B)  - GCA)  IS  -I  - 0.999999998  » -1.999999998 


TYPE  THE  CODE 
?20 

WHAT  IS  A? 

?0 

INPUT  B* 

7 3*  141  592 

HOW  MANY  SUBINTERVALS  ON  CA#BI? 

7200  N 

INTERVAL  t 0 # 3*  141592  I • 200  TRAPEZOIDS 

APPROXIMATION  TO  INTEGRAL  IS  *******  1.999958882  ******* 

GCB)  - GCA)  IS  1 - -1  =2 


TYPE  THE  CODE 
7 20 

WHAT  IS  A? 

70 

INPUT  B- 
710 

HOW  MANY  SUBINTERVALS  ON  CA#BI? 

7200 

INTERVAL  C 0 # 10  ]•  200  TRAPEZOIDS 

APPROXIMATION  TO  INTEGRAL  IS  *******  1.825087846  ******* 

GCB)  - GCA)  1$  0.839071529  - -1  * 1*839071529 


TYPE  THE  CODE 
7 300 


5.4 
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A listing  of  the  program  /FUNDTHEO/. 


1 OPR.  "THIS  PR0GRAM  C0MPUTES  APPR0X IMATI0NS  T0  THE  DEFINITE*’ 

1 2PR • *' INTEGRAL  0F  THE  FUNCTION  Y - F<X>  WHICH  Y0U  SUPPLIED  0N  ” 
1 4PR. "LINE  100#  0VER  AN  INTERVAL  CA#B3*  AND  COMPARES  THE  ” 

1 6PR. "APPR0X IMATI0N  T0  G<B)-G<A>#  WHERE  Z = GCX)  IS  A*1 
1 8 PR • *'PR I M I T I VE  0F  F<X>#  WHICH  YOU  SUPPLIED  0N  LINE  200*** 

20PR. 

22PR. "WHAT  IS  A?”  INPUT  A 
24PR.  ’’INPUT  B.M  INPUT  B 
26  IF  A>-R  G0T0  296 

30  PR.  ”H0W  MANY  SUBINTERVALS  0N  CA#B3?M 

31  INPUT  N 

32  LET  N = INT<ABS<N>> 

40  LET  S=0 

50  LET  Ws<B-A)/N 
60  F0R  X*A  T0  B STEP  W 
100  LET  Y a SI N( X ) 

105  IF  XsB  G0T0  125 
110  IF  X=A  G0T0  i25 
I 15  LET  S=S+C2*Y> 

120  G0T0  130 
125  LET  S = S+Y 
130  NEXT  X 

135  LET  T-<  < B-A) /<  2*N) > *S 

1 40 PR . "INTERVAL  C ’*i  A t ”» **i  B*  **3  • ”1  Nl  "TRAPEZ0 1 DS** 

142  PR. 

1 44PR*”APPR0XIMATI0N  T0  INTEGRAL  IS  ****+**••$  Ti ••***.****•• 

146PR. 

195  LET  L*0 
197  LET  X=A 
200  LET  Z*-C0S<X> 

205  IF  L=!  G0T0  250 

210  LET  L-L+ 1 

215  LET  C=Z 

220  LET  X = B 

225  G0T0  200 

250  LET  D=Z 

255  LET  R-O-C 

260  PR.  **G<  B)  - GCA)  IS**I 01  ••-•*! Cl *’=*’! R 
265  PR. 

270  PR. 

2 75  PR.  "TYPE  THE_C0DE"  INPUT  E 
280  IF  E ■ 20  G0T0  20 

285  IF  E*30  G0T0  30 
290  IF  E=300  G0T0  300 

292  PR. **THE  C0OE  MUST  BE  20>  30#  0R  300. ** 

294  G0T0  275 

296  PR. *’A  MUST  BE  LESS  THAN  B.** 

298  G0T0  22 

300  END  ^ 
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SKETCHING  GRAPHS  OP  FUNCTIONS 


Constructing  the  Graph  of  a Function 

In  this  section  we  are  concerned  with  methods  for 
sketching  graphs  of  functions.  One  approach  to  sketching 
a graph  of  the  function  y — f (x)  is  to  select  various  values 
for  X/  to  calculate  the  corresponding  values  for  y by  sub- 
stituting values  of  x into  the  function  f (x) , and  to  plot  the 
set  of  ordered  pairs  of  numbers  { (x,  f (x) ) } on  rectangular 
coordinate  paper , using  an  appropriate  scale  for  the  x and 
y axes.  When  enough  ordered  pairs  have  been  calculated  and 
plotted  to  enable  you  to  predict  the  shape  of  the  graph  with 
reasonable  accuracy/  the  points  can  be  connected  by  smooth 
lines  (not  necessarily  straight  lines)  to  obtain  a graph  of 
the  function  y = f (x) . If  the  graph  has  been  precisely  con- 
structed/ it  is  possible  to  obtain  much  information  about 
the  function  by  studying  the  graph. 

When  analysing  a function/  one  should  consider  the  be- 
havior of  the  function  relative  to  any  open  or  closed  interval 
in  the  domain  of  the  function  with  respect  to  the  following 
mathematical  concepts: 

(1)  Is  the  function  increasing  or  decreasing? 

(2)  Is  the  function  one-to-one? 

(3)  Is  the  function  onto? 

(4)  Does  the  function  have  an  inverse  function? 

(5)  Does  the  function  have  any  real  zeros  (roots)? 

(6)  Is  the  function  strictly  increasing  or  strictly 
decreasing? 

(7)  Is  the  function  constant? 

(8)  Does  the  function  have  maximum  or  minimum  values? 

(9)  Does  the  function  have  relative  maximum  or  relative 
minimum,  values? 

(10)  Is  the  function  concave  upward  or  concave  downward? 

(11)  Are  there  any  real  numbers  for  which  the  function 
is  undefined? 

(12)  Are  there  any  real  nurjbers  for  which  the  function 
is  discontinuous? 
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(13)  Are  there  any  real  numbers  for  which  the  function 
does  not'  have  a derivative? 

(14)  What  is  the  range  of  the  function? 

(15)  What  is  the  domain  of  the  function? 

Check  your  sources  of  information  (books,  other  students, 
teachers,  etc,)  to  make  certain  that  you  understand  each  math- 
ematical term  that  is  underlined. 


Problem  24;  Write  computer  programs  which  will  accept  as  input 
numerical  values  of  x and  give  as  output  values  of  f (x)  for 
each  of  the  following  functions . Use  your  program  to  obtain 
enough  points  on  the  graph  of  each  function  so  that  you  can 
plot  an  accurate  graph  of  the  function. 


A,  f (x)  = x2-3x+2 

B,  f (x)  = -x2+3x-2 

C,  f (x)  = cos  (x) 

D,  f (x)  * 2cos(x) 

E,  f (x)  = cos(2x) 

P,  f (x)  = log  x 

H . 

G.  f (x)  = Jsin  (i)  if  x*0 

L 0 if  x=0 

H.  f (x)  = x(sin(x)) 

ex 

I.  f (X)  = S_ 


38 
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This  module  is  a guessing  game*  The  program  called 
* /GRAPHl/  contains  five  different  functions  of  x*  For  any 
of  the  functions  and  for  any  real  value  of  x which  you 
select/  the  program  will  give  as  output: 

(1)  that  value  of  x, 

(2)  the  value  of  the  function/ 

(3)  the  value  of  the  derivative  of  the  function/ 

(4)  the  value  of  the  second  derivative  of  the  function* 

Your  job  is  to  select  as  many  of  these  functions  as 
you  want,  sketch  graphs  of  the  functions  and  their  first 
and  second  derivative  functions/  and  try  to  guess  what  the 
functions  are  and  what  the  derivative  functions  are.  To 
sketch  the  graphs  you  will  have  to  make  shrewd  selections 
for  your  values  of  x in  order  to  get  points  which  will  assist 
you  in  obtaining  accurate  graphs. 

You  should  try  to  find  maximum  and  minimum  points  on  the 
graphs,  points  where  the  functions  are  discontinuous  or  unde- 
fined/ intervals  where  the  functions  are  either  increasing  or 
decreasing,  intervals  where  the  functions  are  either  concave 
upward  or  downward  and  points  where  the  graphs  cross  the  x 
and  y axes. 

After  you  have  loaded  the  program  /GRAPHl/  and  typed: 
>RUN 

you  will  get  the  message 

NOW  TYPE  THE  CODE  NUMBER  OF  THE  FUNCTION 
THAT  YOU  WANT  TO  WORK  WITH. 

* 

The  code  numbers  of  the  functions  are  100,  200,  300, 

400,  and  500;  so  you  can  select  a function  by  typing  any 
one  of  these  numbers  following  the  question  mark.  After 
you  have  selected  a function  you  will  be  asked  to  type  a 
vaxue  of  x.  When  you  supply  the  value  of  x you  will  ob- 

O 
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tain  as  output  that  value  of  x,  together  with  the  values 
of  the  function  and  its  first  and  second  derivatives*  The 
message : 

TYPE  A VALUE  FOR  X 

will  be  repeated  five  times,  then  you  will  have  the  option 
of  getting  values  of  the  function  and  its  derivatives  for 
five  more  values  of  x or  of  going  on  to  another  function* 

A sample  output  from  this  program  is  shown  on  the  next 

page. 
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A sample  run  of  the  program  /GRAPHl/. 


MOW  TYPE  THE 

CODE 

[ NUMBER  0F 

THE 

FUNCTION 

THAT 

YOU  WANT  TO 

WORK 

WITH. 

? 100 

TYPE 

A VALUE 

FOR 

X. 

?-2 

X = 

-2 

Y = 

4 

01 

s 

-4 

02  = 

2 

TYPE 

A VALUE 

FOR 

X. 

? -1 

X * 

-1 

Y * 

1 

01 

-2 

02  = 

2 

TYPE 

A VALUE 

FOR 

X. 

?0 

X = 

0 

Y = 

0 

01 

s 

0 

02  * 

2 

TYPE 

A VALUE 

FOR 

X. 

? 1 

X = 

1 

Y = 

1 

01 

s 

2 

02  = 

2 

TYPE 

A VALUE 

FOR 

X. 

# 

2 

X = 

2 

Y » 

4 

01 

s 

4 

02  = 

2 

IF  Y0U  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION* 
TYPE  ST0PJ  FOR  MORE  INFORMATION*  TYPE  60. 

? STOP 

TO  OBTAIN  ACCESS  T0  ANOTHER  FUNCTION  TYPE  ITS  CODE 
NUMBER*  100*  200*  300*  400*  OR  500*  TO  STOP  THE 
PROGRAM*  HIT  THE  ESCAPE  KEY. 

NOW  TYPE  THE  CODE  NUMBER  OF  THE  FUNCTION 


THAT  YOU  WANT  TO 

WORK 

WITH. 

? 400 

TYPE  A VALUE 

FOR 

X* 

? -5 

X = -5 

Y * 

0.9589C  4275 

01  = 

0.283662185 

02  = 

-0.958924275 

* 

TYPE  A VALUE 

o 

F0R 

X* 

t 

-2 

X * -2 

Y « 

-0*909297427 

01  s 

-0.416146837 

02  = 

0.909297427 

TYPE  A VALUE 

F0R 

X* 

?0 

X * 0 

Y * 

0 01  * 

1 

02  s 0 

TYPE  A VALUE 

F0R 

X* 

? 1 

- 

■ --  - 

- - - 

X * 1 

Y * 

0.841470985 

01  = 

0.540302306 

02  = 

-0*841470985 

TYPE  A VALUE 

F0R 

X* 

? 10 

X = 10 

Y * 

-0*544021111 

01  S 

-0.839071529 

D2  = 

0*544021111 

IF  YOU  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION* 

TYPE  ST0PJ  IF  YOU  WANT  MORE  INFORMATION*  TYPE  G0* 

? STOP 

TO  OBTAIN  ACCESS  T0  ANOTHER  FUNCTION  TYPE  ITS  CODE 
NUMBER*  100*  200*  300*  400*  OR  500.  TO  STOP  THE 
PROGRAM*  HIT  THE  ESCAPE  KEY. 

NOW  TYPE  THE  CODE  NUMBER  0F  THE  FUNCTION 
THAT  YOU  WANT  T0  WORK  WITH. 

? 300 

TYPE  A VALUE  FOR  X. 

? -4 

X = -4  Y * 0.018315639  01  » 0.018315639  02  * 

0.018315639  A 

TYPE  A VALUE  FOR  X.  91 

?*-**ESC*  305 


A listing  of  the  program  /GRAPHl/.  g ^ 

3 4 PRINT  “ “ 

36  PRINT  “NOW  TYPE  THE  C0DE  NUMBER  0F  THE  FUNCTION** 

38  PRINT  “THAT  Y0U  WANT  T0  W0RK  WITH**" 

44  INPUT  C 

46  IF  C = 100  60T0  100 
48  IF  C s 200  60T0  200 
50  IF  C - 300  60T0  300 
52  IF  C - 400  60T0  400 
54  IF  C - 500  60T0  500 

56  PRINT  “YOU  MUST  TYPE  100,  200,  300,  400,  0R  500. ** 

58  PRINT  “TRY  AGAIN. “ 

60  G0T0  44 
100  LET  N = 0 

102  PRINT  “TYPE  A VALUE  F0R  X*“ 

105  INPUT  X 
107  LET  Y=X*X 
1 10  LET  D1=2*X 
1 15  LET  02=2 

120  PRINT  “X  = “1X1“  Y = “1Y1“D1  = “!D!!“  D2  * “1D2 
1 25  LET  N = N+l 
130  IF  N=5  G0T0  150 
135  G0T0  102 

150  PRINT  “IF  Y0U  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION,** 
155  PRINT  “TYPE  STOP!  FOR  MORE  INFORMATION,  TYPE  G0*“ 

165  INPUT  T$ 

170  IF  IS  <T$,"G0“,O>  GOTO  100 
175  IF  IS  < TS, “STOP**, 0>  GOTO  600 

180  PRINT  “YOU  MUST  TYPE  STOP  OR  GO*  TRY  AGAIN. “ 

185  GOTO  165 
200  LET  N=0 

202  PRINT  “TYPE  A VALUE  FOR  X*“ 

205  INPUT  X 
207  LET  Y=X*<X-4> 

210  LET  D!=2*<X-2> 

215  LET  D2=2 

220  PRINT  “X  = **IXI“  Y * “|Y|“  D!  = “|D!|“  D2  = “|D2 

225  LET  N=N+1 

230  IF  N * 5 GOTO  250 

235  GOTO  202 

250  PRINT  “IF  YOU  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION,** 
255  PRINT  “TYPE  STOP!  IF  YOU  WANT  MORE  INFORMATION,  TYPE  GO*** 

260  INPUT  TS  . _ - - — 

270  TF  IS  <T$,“G0“,O>  GOTO  200 
275  IF  IS  <T$,“ST0P”,O>  GOTO  600 

280  PRINT  “ YOU  MUST  TYPE  STOP  0R  60*  TRY  AGAIN. “ 

285  GOTO  260 
300  LET  N = 0 

302  PRINT  “TYPE  A VALUE  FOR  X*“ 

305  INPUT  X 

307  LET  Y * EXP(X) 

310  LET  D 1=Y 
3 1 5 LET  D2=Y 

320  PRINT  “X  a “1X1“  Y = “|Yl“  01  * “iDll**  D2  * “|D2 

325  LET  N = N+l 

330  IF  N = 5 GOTO  350 

335  GOTO  302 

350  PRINT  **  IF  YOU  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION, 
355  PRINT  “TYPE  STOP!  IF  YOU  WANT  MORE  INFORMATION,  tYPE  G0.“ 
360  INPUT  TS 

370  IF  IS  (TS,“G0“,O>  GOTO  300 
375  IF  IS  ( TS , “STOP**, 0)  GOTO  600 

380  PRINT  “ YOU  MUST  TYPE  STOP  OR  GO.  TRY  AGAIN*** 

385  GOTO  360  --  je-  • . 


A listing  of  the  program  /GRAPHl/  (continued) 


400  LET  N=0 

402  PRINT  "TY PE  A VALUE  FOR  X." 

405  INPUT  X 

407  LET  Y = SIN(X) 

410  LET  D1*C0S(X) 

415  LET  02  *-SIN<X> 

420  PRINT  **X  * "IXI"  Y * "IYI"  01  '*  "IOU"  02  = "102 

425  LET  N * N+l 

430  IF  N * 5 60 TO  450 

435  GOTO  402 

450  PRINT  **IF  YOU  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION 
455  PRINT  "TYPE  ST0P1  IF  YOU  WANT  MO RE  INFORMATION#  TYPE  GO*" 
460  INPUT  TS  ' v 

470  IF  IS  <TS#"G0"#O)  GOTO  400 
475  IF  IS. <TS#"ST0P"#O>  GOTO  600 

480  PRINT  " YOU  MUST  TYPE  STOP  OR  G0.  TRY  AGAIN." 

485  GOTO  460 
500  LET  N=0 

502  PRINT  "TYPE  A VALUE  FOR  X." 

505  INPUT  X 

507  LET  Y = X*(X*X-X+1) 

510  LET  01-X*<  3*X-2)+l 
515  LET  02=6*X~2 

520  PRINT  "X  = "IX*"  Y = "IYI"  01  = "IOU"  02  * "102 

525  LET  N*N+1 

530  IF  N = 5 GOTO  550 

535  GOTO  502 

550  PRINT  "IF  YOU  HAVE  ENOUGH  INFORMATION  ABOUT  THIS  FUNCTION 
555  PRINT  "TYPE  STOP*  IF  YOU  WANT  MORE  INFORMATION# TYPE  GO." 
560  INPUT  TS  _ 

570  IF  IS  (TS#"G0"#O)  GOTO  500 
575  IF  I S< TS# •#ST0 P"# 0)  GOTO  600 

580  PRINT  " YOU  MUST  TYPE  STOP  0R  GO.  TRY  AGAIN*" 

585  GOTO  560 

600  PRINT  "TO  OBTAIN  ACCESS  TO  ANOTHER  FUNCTION  TYPE  ITS  C0OE 
605  PRINT  "NUMBER#  100#  200#  300#  400#  0R  500.  TO  STOP  THE" 
610  PRINT  "PROGRAM#  HIT  THE  ESCAPE  KEY." 

650  GOTO  36 


This  module  called  /GRAPH2/  is  a modification  of  /GRAPH1/. 

It  is  a game  for  guessing  functions;  however  your  guesses  must 
be  based  upon  less  information  than  is  given  as  output  from 
/GRAPHl/.  The  program  /GRAPH2/.  contains  the  first  and  second 
derivative  functions  for  six  different  functions  of  x.  For 
any  of  these  six  functions  and  for  any  real  value  of  x which 
you  select/  the  program  will  give  you  for  output: 

(1)  the  value  of  x which  you  selected/ 

(2)  the  value  of  the  derivative  of  the  function/ 

(3)  the  value  of  the  second  derivative  of  the  function. 

Your  problem  is  to  select  some  or  all  of  the  six  functions/ 
to  get  as  much  information  about  the  first  and  second  derivatives 
of  the  functions  as  you  need  in  order  to  sketch  accurate  graphs 
of  the  derivative  functions/  and  finally  to  guess  the  functions 
whose  first  and  second  derivative  functions  are  stored  in  the 
computer . 

You  should  review  the  information  concerning  the  relation- 
ship of  a function  to  its  derivative  function/  and  be  especially 
concerned  with  maximum  and  minimum  points/  inflection  points / 
critical  points,  and  intervals  where  the  function  is  increas- 
ing or  decreasing  and  concave  upward  or  downward. 

After  you  have  loaded  /GRAPH 2/  and  typed: 

>RUN 

you  will  get  the  message 

NOW  SELECT  A FUNCTION  BY  TYPING  ITS  CODE 
NUMBER;  100,  200,  300,  400,  500,  OR  600. 

9 

* 

You  can  now  select  any  of  the  six  functions  by  typing 
its  code  number  following  the  question  mark.  After  you  have 
chosen  a function  to  study,  you  will  be  asked  to  type  a value 
for  x.  After  typing  a value  for  x you  will  be  given  as  output 
that  value  of  x together  with  corresponding  values  of  the  first 
and  second  derivative  functions  of  the  function  which  you  are 
attempting  to  guess. 


* to 


* * 
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The  message 

TYPE  THE  VALUE  FOR  X. 

will  be  repeated  five  times#  then  you  will  have  the  option 
of  getting  values  of  tihe  derivatives  for  five  more  values 
of  x or  of  going  on  to  another  function.  When  you  have 
enough  information  about  the  functions#  hit  the  escape  key* 

A sample  output  from  this  program  is  shown  on  the  next 

page. 
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A sample  output  of  the  program  /GRAPH 2/* 
N0W  SELECT  A FUNCTION  BY  TYPING  ITS  CODE 


6*10 


••  *# 


NUMBER! 

100*  200* 

300* 

400*  5A0*  OR 

600* 

7 200 

TYPE 

THE 

VALUE  OF 

X* 

7-1 

X* 

-1 

FI* 

0*367879441 

F2* 

0*367879441 

TYPE 

THE 

VALUE  OF 

X* 

70 

X* 

0 

FI* 

1 F2* 

1 

. TYPE 

THE 

VALUE  OF 

X* 

72 

X* 

2 

FI* 

7*389056099 

. F2* 

7*389056099 

TYPE 

THE 

VALUE  OF 

X* 

73 

- 

i X* 

3 

FI* 

20*08553692 

F2* 

20*08553692 

! TYPE 

THE 

VALUE  OF 

X- 

i ? -2 

t x= 

-2 

FI* 

0*135335283 

F2* 

0*135335283 

1 0 ASSIGN  5 MORE  VALUES  TO  X*  TYPE 

GO* 

4*48168907 


7*389056099 


2980*957987 


F2* 


F2* 


F2= 


TO  LEAVE  THIS  FUNCTION*  TYPE  STOP* 

7 GO 

TYPE  THE  VALUE  OF  X* 

7 1*5 

X*  1*5  FI* 

TYPE  THE  VALUE  OF  X* 

? 2 

X = 2 FI* 

TYPE  THE  VALUE  OF  X* 

78 

X*  8 FI* 

TYPE  THE  VALUE  OF  X* 

7-9 

X*  -9  FI*  1 * 23  4098041 E -04  F2= 

TYPE  THE  VALUE  OF  X* 

77 

X*  7 FI* 

TO  ASSIGN  5 MORE  VALUES  TO  X*  TYPE  GO* 

TO  LEAVE  THIS  FUNCTION*  TYPE  STOP* 

?STP 

TYPE  EITHER  STOP  OR  GO*  TRY  AGAIN* 

7 STOP  - 

TO  GET  ANOTHER  FUNCTION  TYPE  ITS  CODE  NUMBER* 
TO  STOP  THE  PROGRAM  HIT  THE  ESCAPE  KEY* 

7 300 


1 * 23 4098041 E -04 

1096*633158  F2* 


4*48168907 

7*389056099 

2980*957987 

1 *23409804 1E~04 

1096*633158 


TYPE 

A 

VALUE 

FOR 

X* 

7-1 

■ 

X*  , 

\ 

FI  * 

1 

F2* 

tyre 

A 

VALUE 

FOR 

X* 

70  ' 

X*0* 

FI 

AND  1 

F2  00 

i NOT 

EXIST 

TYPE 

A 

VALUE 

FOR 

X* 

73 

X* 

3 

FI* 

1 

F2= 

TYPE 

A 

vALUE 

FOR 

X* 

75 

X * 

5 

FI* 

1 

F2* 

TYPE 

A 

VALUE 

FOR 

X* 

78 

X* 

8 

FI* 

1 

F2* 

TO  ASSIGN  5 

MORE 

VALUES  TO 

X* 

TO  LEAVE  THIS  FUNCTION*  TYPE  STOP* 
? **-**ESC?  510 


i \ 

,46 


* 


I I 


*1/ 


LC 
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A listing  of  the  program  /GRAPH2/. 


46 
*n, 
50 
52 
54 
56 
58 
60 
62 
64 
66 
68 
70 
1 00 
1 05 
t to 
1 t 5 
t 20 
1 25 
1 30 
t 35 
1 40 
t 50 
t 55 
t 60 
t 65 
t 7 0 

1 75 
t 80 
200 
205 
210 
215 
220 
225 

2 30 
235 
2 40 
2 50 
255 
2 60 
2 65 

2 70 
275 
280 

3 00 
3 05 
3 10 
3 15 
320 
325 
3 30 
335 
340 
3 50 
355 


••  «• 


PR. 

PR. 

PR. 

INPUT  C 
IF  C = 100 


"NOW  SELECT  A 
"NUMBER!  100* 


FUNCTION 
200*  300* 


BY  TYPING  ITS  CODE 
400*  500*  OR  600* 


•• 


IF 

IF 

IF 

IF 

IF 

PR. 

PR. 


= 201 
= 300 
* 400 
- 500 
a 600 
"YOU 
"500* 


GOTO  100 
GOTO  200 
GOTO  300 
GOTO  400 
GOTO  500 
GOTO  600 
MUST  TYPE 
OR  600* 


EITHER  100* 
TRY  AGAIN*" 


200*  300*  400* 


GOTO  52 
LET  N = 0 

PR.  "NOW  TYPE  A VALUE  FOR  X." 

INPUT  X 
LET  F 1 *2 
LET  F2=0 

PR.  "X  * "IX!"  FI  * "!F1I"  F2  * ”|F2 
LET  N * N+l 
IF  N=5  GOTO  150 
GOTO  105 

PR*  "TO  ASSIGN  5 MORE  VALUES  TO  X*  TYPE 
PR.  "TO  GET  ANOTHER  FUNCTION  OF  X*  TYPE 
INPUT  TS 


GO*" 

STOP- 


IS  <T$*"G0"*O>  GOTO 
IS  < T$* "STOP"* 0>  GOTO 
"YOU  MUST  TYPE  GO  OR 


too 

700 
STOP. 


F 1*  -f"l  F 1 1 " F2*  "IF2 


250 


IF 
IF 
PR* 

GOTO  160 
LET  N*0 

PR.  "TYPE  THE  VALUE  OF  X. 
INPUT  X 

let  fi*e::p<X) 

LET  F2*F! 

PR.  *‘X=  "IX!" 

LET  N»N+! 

IF  N*5  GOTO 
GOTO  205 

PR*  "TO  ASSIGN  5 MORE  VALUES 
PR*  "70  LEAVE  THIS  FUNCTION* 
INPUT  TS 

IF  IS  <T$*"G3"*0>  GOTO  200 
IF  IS  < TS* "STOP"* 0)  GOTO  700 
PR.  "TYPE  EITHER  STOP  OR  G0. 
GOTO  260 
LET  N=0 

PR.  "TYPE  A VALUE  FOR  X." 
INPUT  X 

IF  X*0  GOTO  385 
IF  X<0  LET  Fl=- 1 
IF  X' 0 LET  F 1* 1 
LET  F iisO 


TRY  AGAIN. 


•» 


TO  X*  TYPE  GO* 
TYPE  STOP*" 


TRY  AGAIN* 


•t 


PR.  "X*  "IX!" 
LET  N » N+l 
IF  N* 5 GOTO  360 
GOTO  305 


Fl=  "Jrl!"F2s 


"JF2 


A listing  of  the  program  /GRAPH 2/  (continued) 

^ * * 

360  PR*  MT0  ASSIGN  5 MORE  VALUES  70  X#  TYPE  GO." 

365  PR*  M T0  LEAVE  THIS  FUNCTION#  TYPE  STOP.** 

370  INPUT  T$ 

372  IF  IS  <T$#"G0’*#O>  G0T0  300 

374  IF  IS  ( T$#  "STOP**#  0)  G0T0  700 

376  PR*  MTYPE  EITHER  G0  0R  STOP*  TRY  AGAIN*** 

380  GOTO  370 

385  PR.  **X=0#  FI  AND  F2  00  NOT  fe^-ST." 

390  GOTO  340 
400  LET  M=0 

405  PR.  **TYPE  A VALUE  FOR  X*M 

410  INPUT  X 

415  LET  F!*-SIN(X> 

420  LET  F2=-C0S(X) 

425  PR.  "X*  "IX!"  FI*  *'l F 1 1 "F2®  "I  F2 

430  LET  N=N+I 

435  IF  N*5  GOTO  450 

440  GOTO  405 

450  PR*  "TO  GIVE  X FIVE  MORE  VALUE.S#  TYPE  GO." 

455  PR.  "TO  GET  ANOTHER  FUNCTION#  TYPE  STOP." 

460  INPUT  T$ 

465  IF  IS  <T$#"G0"#O>  GOTO  400 

470  IF  IS  <TS#"ST0P"#0>  G0T0  700 

475  PR.  " TYPE  STOP  OR  GO*  TRY  AGAIN*** 

480  GOTO  460 
500  LET  N=0 

505  PR.  "TYPE  A VALUE  FOR  X." 

510  INPUT  X 
515  LET  FI=3*X*X 
520  LET  F2*6*X 

525  PR.  "X*  "I  XI"  FI=  "I F 1 1 "F2*  "I F2 
530  LET  N=M+! 

535  IF  N*5  GOTO  550 
540  GOTO  505 

550  PR.  "TO  GIVE  X MORE  VALUES#  TYPE  G0.*' 

555  PR.  " TO  GO  TO  A DIFFERENT  FUNCTION#  TYPE  STOP." 
560  INPUT  TS 

565  IF  IS  ( TS# "GO**# 0)  GOTO  500 

570  IF  IS  (TS# "STOP**'#  0)  G0T0  700 

575  PR.  "TYPE  EITHER  STOP  OR  GO.  TRY  AGAIN.** 

580  GOTO  S60 
600  LET  N*0 

605  PR.  "TYPE  A VALUE  FOR  X*" 

'610  INPUT  X 

615  LET  F 1 *X*( X-3 ) *2 

620  LET  F2=2*X-3 

625  PR.  "X*  "IX!"  FI*  "!"1!**F2  * **!  F2 

630  LET  N*N+I 

635  IF  N=5  GOTO  650 

64C  GOTO  605 

650  PR.  "TO  GIVE  MORE  VALUEf  TO  X#  TYPE  GO." 

655  PR.  "TO  G0  TO  ANOTHER  FUNCTION#  TYPE  STOP.** 

660  INPUT  T$ 

665  IF  IS  ( TS# "GO**# 0)  GOTO  600 

670  IF  IS  (T$#**ST0P**#0)  G0T0  700 

675  PR.  "YOU  MUST  TYPE  GO  0R  STOP.  TRY  AGAIN." 

680  GOTO  660 

700  PR.  "TO  GET  ANOTHER  FUNCTION  TYPE  ITS  CODE  NUMBER. 
705  PR.  "TO  STOP  THE  PROGRAM  HIT  THE  ESCAPE  KEY." 

710  GOTO  52  *0 


